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The behavior of polyelectrolytes between charged surfaces immersed in semi-dilute solutions 
is investigated theoretically. A continuum mean field approach is used for calculating numerically 
concentration profiles between two electrodes held at a constant potential. A generalized contact 
theorem relates the inter-surface forces to the concentration profiles. The numerical results show 
that over-compensation of the surface charges by adsorbing polyelectrolytes can lead to effective 
attraction between equally charged surfaces. Simple scaling arguments enable us to characterize 
qualitatively the inter-surface interactions as function of the fraction of charged monomers p and 
the salt concentration q,. In the low salt regime we find strong repulsion at short distances, where 
the polymers are depleted from the inter-surface gap, followed by strong attraction when the two 
adsorbed layers overlap. The magnitude of this attraction scales as p 1 ^ 2 and its dominant length 
scale is proportional to a/p 1 ^ 2 , where a is the monomer size. At larger distances the two adsorbing 
surfaces interact via a weak electrostatic repulsion. For strong polyelectrolytes at high salt concen- 
tration the polymer contribution to attraction at short distances scales as p/c^ 2 and the length 
scale is proportional to tt s a 2 /p, where k^ 1 is the Debye-Hiickel screening length. For weak polyelec- 
trolytes at high salt concentration the interaction is repulsive for all surface separations and decays 
exponentially with a decay length equal to kJ 1 . The effect of irreversible adsorption is discussed as 
well and it is shown that inter-surface attraction can be obtained in this case as well. 
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I. INTRODUCTION 

Polymers are known to affect the interactions of colloidal particles in solution . Adsorption of charged poly- 
mers {polyelectrolytes) to oppositely charged colloids may turn inter-colloidal repulsion into attraction, leading to 
fiocculation. This phenomenon is used in industrial applications such as water filtration, paper making and mineral 
processing. The reversed process is useful as well, since adsorbed polyelectrolytes (in different conditions) can also 
stabilize colloidal suspensions such as paint, ink or medical suspensions against attractive forces (e.g., van der Waals 
forces) . 

One of the most common techniques to study experimentally the adsorption of polyelectrolytes between two surfaces 
is the Surface Force Apparatus (SFA) which allows delicate measurements of inter-surface forces at distances as 
small as a few Angstroms. In these experiments [p|-^2[, attractive and repulsive forces have been observed, depending 
on specific details, such as the type of polyelectrolyte, its concentration, the ionic strength of the solution, etc. In other 
experiments it is possible to measure the disjoining pressure of thin liquid films in the presence of polyelectrolytes 
and as a function of the film width Jl3[ . Both repulsive and attractive forces have been measured using this method. 

Adsorption of polyelectrolytes was treated theoretically in a discrete model |l4|-|l7t , where the state of the system 
is described by occupation fractions of monomers, ions and solvent molecules on a discrete lattice. Within mean field, 
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the equilibrium state of the system corresponds to the maximal contribution to the partition function of the system, 
and can be calculated numerically. Bohmer et al. [Q have used this model to calculate force curves at relatively 
short distances (up to 30 molecular layers). In addition, Monte Carlo computer simulations of polyelectrolytes 
between flat surfaces ]l9[ ] and between charged spheres ppf ] provide valuable hints concerning the complex behavior 
of polyelectrolytes. However, they are limited to relatively short chains and small inter-surface distances due to 
computation time limitations. 



Another theoretical approach is a continuum one 2 1 - 2q] where the concentration of the different species are 



taken to be continuous functions of the spatial coordinates. The mean-field state can be calculated by solving 
two differential equations for the polymer concentration and the electrostatic potential derived using a variational 
procedure. Varoqui et al. J22|,^3| used the continuum approach to investigate polyelectrolyte adsorption onto one 
surface, while Podgornik [ p7| used a similar formalism to calculate inter-surface forces. In those works the non-linear 
excluded volume interaction between the monomers has not been considered. Chatellier and Joanny |p8| used a 
linearized version of a similar approach to study the inter-surface interactions for polyelectrolytes in a poor solvent. 
Recently, we have been able to derive some simple scaling relations for the adsorption of polyelectrolytes onto a 
single charged surface |?6| . These scaling relations were compared to the exact numerical solutions of the differential 
equations and to existing experimental results. The agreement was reasonable in two opposite limits: (i) low salt 
concentration (no electrostatic screening) and (ii) high salt concentration (strong screening). 

In the present work, the continuum model is used to study polyelectrolyte adsorption between two parallel surfaces. 
The advantage of our model is that the connectivity of the polymer chains, the excluded volume repulsion between 
monomers in a good solvent and the Coulomb interactions between the charged monomers, counter-ions, co-ions 
and surface charges are all taken into account. The mean field equations are solved numerically in order to obtain 
concentration profiles, and the inter-surface interactions are then calculated from the free energy. In addition we 
extend our earlier scaling approach of one adsorbing surface |26| ] to the case of two interacting surfaces. The result is 
a qualitative description of the inter-surface interactions as function of the polyelectrolyte charge and the amount of 
salt in the solution. 

The paper is organized as follows: in the next section we describe the mean-field approach. In Sec. Ill we present 
numerical results obtained from solving the mean field equations, and in Sec. IV we use simple scaling arguments to 
describe the inter-surface interactions. In Sec. V we study the effect of irreversible adsorption both numerically and 
analytically and in Sec. VI we compare our results with experiments. Finally we present our conclusions and some 
future prospects. 

II. THE MEAN FIELD APPROACH 
A. The Basic Equations 

The model system consists of a semi-dilute solution of polyelectrolytes in a good solvent placed between two flat 
surfaces (Fig. 1). The solution contains charged polymer chains, counter-ions and a monovalent electrolyte (salt). 
Having in mind the experimental setup of the surface force apparatus (SFA) discussed below, we consider a system 
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which is coupled to a bulk reservoir of polyelectrolyte chains and salt. 

As discussed elsewhere |24j,^5|, the charge distribution along the chains depends on the type of polyelectrolyte as 
well as on the local conditions, such as the pH and the electrostatic potential. However, at low electrostatic potentials 
\/3eip\ <C 1, where ip is the electrostatic potential, [3 = 1/ksT the inverse thermal energy and e the electron charge, 
the differences between the various charge distribution models are small. We therefore assume hereafter a uniform 
charge distribution along the polymer chains with a fractional charge pe attached to each monomer. 

In the mean field approach, the free energy of the system is expressed in terms of the local electrostatic potential 
V'(r) at a point r and the polymer order parameter <f>(r) which is related to the local monomer concentration through 
Cm(r) = |</>(r)| 2 . The relation between the polymer order parameter and the monomer concentration is analogous to 
the relation between the wave function and the probability density of a particle in quantum mechanics. The excess 
free energy with respect to the bulk can be divided into three contributions P2|-f25|: 



F = I /(r)dr = / { / pol (r) + / ions (r) + / el (r) f-dr (1) 

The polymer contribution is 



/poi(r) = k B T 



i-|V0| 2 + ±t# 4 -#) 



ihiF-K) (2) 



where the first term is the polymer response to local variations of the concentration and is due to the connectivity 
of the polymer chain, a being the effective monomer size. The second term represents the short ranged monomer- 
monomer interaction and can be viewed as representing an effective volume of a single monomer. For a polymer in 
a good solvent v is positive. However, since v represents an effective monomer-monomer interaction it can also be 
negative (in a poor solvent) or zero (in a theta solvent) requiring higher order terms in 4> 2 to be included in the free 
energy. For simplicity we will limit ourselves to good solvent conditions but the formalism can be easily generalized 
to other conditions as well. The last term couples the system to a reservoir, fi p being the polymer chemical potential 
and 4>\ the bulk monomer concentration. 

The non-electrostatic contribution of the small (monovalent) ions is due to their translation entropy and is equal 

to 



/ioBsW = Y,{ kBT c * lncl ~ <M4M - 4) - - 4)} 



(3) 



i=± 

where c l (r) is the local concentration of the i = ± ions (cations and anions) and cl, [i % are the bulk concentration 
and chemical potential, respectively. In the most general case the solution contains two types of negative ions: the 
counter- ions which dissociate from the polymer chains and the salt anions. In the reservoir the concentration of 
negative ions has two contributions = Cb +p4>l where c& is the electrolyte bulk concentration, while for the positive 
ions cjj~ = q,. In principle, one could consider the two types of negative ions separately, but for clarity we take the 
two types of negative ions to be identical. 
Finally, the electrostatic contribution is 

/ei(r) = pe^ty + ec+^ - ec'^j - ^-|VVf (4) 

07T 
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The first term is the electrostatic energy of charged monomers. The next two terms represent the positive and negative 
ions, respectively, and the last term is the self energy of the electric field where e is the dielectric constant of the 
solution. The sum of the electrostatic contributions can be integrated by parts using the Poisson equation (derived 
below) and yields F c i = (e/8n) J |V?/>| 2 dr , as expected, plus electrostatic surface terms. 

Minimization of the free energy with respect to c^, (f> and ip yields a Boltzmann distribution for the concentration 
of the small ions, c ± (r) = exp(=F/3e?/>), and two coupled differential equations for cf> and ip p5[ : 

V 2 V(r) = — c b smh((3etP) - — ( P 4> 2 - jxfieP**) (5) 
e e 

^ V 2 0(r) = v{cf - fy) + P 4>(3e^ (6) 
6 

Equation is a generalization of the Poisson-Boltzmann equation including the free ions as well as the charged 
polymers. The first term represents the salt contribution and the second term is due to the charged monomers and 
their counter-ions. Equation |]is a generalization of the self-consistent field equation of neutral polymers p9j . In the 
bulk the potential and the polymer concentration have constant bulk values given by ip — and <j> = cfn,, as can be 
seen in the above equations. 



B. Two Interacting Surfaces 



The interaction of two charged surfaces in a solution containing only small ions (electrolyte) without charged poly- 
mers is well established within the framework of the Poisson-Boltzmann equation f3(i[| . The electrostatic interaction 
between two identically charged surfaces is found to be repulsive within this mean-field like theory However, 
the addition of polyelectrolytes to the solution changes the picture in a subtle way. Experiments |p|~|l3|] show that 
polyelectrolytes reduce this repulsion and might even cause mutual attraction between the two surfaces. 

For simplicity, the surfaces are taken as flat, homogeneous and parallel to each other in order that the physical 
quantities will depend only on the position x between the surfaces (see Fig. 1). The effect of the surfaces is introduced 
through the boundary conditions on the polymer order parameter <p{x) and the electrostatic potential ip(x). In this 
work, both surfaces are assumed to be kept at the same constant potential: 

Ms = A (7) 

and no monomers are adsorbed on the surfaces 

4>\ s = (8) 

Other boundary conditions could have been considered as well. For example, if a fixed surface charge a is assumed 
then the electrostatic boundary condition would include the electric field: tp'\ s = —Ana/e. In real systems neither 
the surface potential nor the surface charge are fixed. The choice of one or another is only an approximation whose 
quality depends on the details of the experimental system. Similarly, for the polymer boundary conditions one could 
consider an adsorbing surface instead of a non-adsorbing one [ ^3| . 
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Given these boundary conditions, the Poisson-Boltzmann and self-consistent field equations uniquely deter- 
mine ip(x) and <f>(x). However, experiments usually probe global properties such as the amount of monomers adsorbed 
per unit area or the inter-surface interactions. 

The total amount of monomers (per unit area) between the two surfaces T(w) can be easily calculated from the 
polymer concentration profile since 

/w/2 
4> 2 {x)dx (9) 
-w/2 

Another measure for the strength of the adsorption is the average monomer concentration divided by the bulk 
concentration: 



i r' 2 <p{x) A _v{w) 



" / - / xl dx = — IT ( 10 ) 

The latter quantity relates to the strength of the adsorption only at small distances. As demonstrated by the numerical 
examples below, at larger distances, V saturates to a constant and the average concentration decreases as 1/w. 

In addition, it is of interest to calculate the total amount of charge (per unit area) carried by the adsorbed polymers 
<j p (w) = peT(w) as compared to the induced surface charge density a s (w) (on a single surface). The latter depends 
on the inter-surface distance as we have chosen to work with constant surface potentials rather than constant surface 
charges. 

The adsorption of polyelectrolytes strongly affects the inter-surface interactions. The excess free energy per unit 
area for two surfaces at a distance w apart can be calculated from the concentration profiles <f>(x) and ij)(x): 

AF(w) = / f[<l>(x),il>(x)]dx -2F 1 (11) 

J -w/2 

where fix) was introduced in eqs. and F{ w — » oo ) = 2F\ is the free energy of two isolated surfaces at infinite 
separation. 

The variation of this free energy with respect to the inter-surface distance w gives the inter-surface pressure (or 
force per unit area): 

It can be shown from eqs. p"l| , |l2| and from 6(AF)/S<p(x) = S(AF) / '5ip(x) = that 

IL(w) = -f{x = 0) (13) 

where fix = 0) is the free energy density (per unit volume) at the mid-plane. This relation is a generalization of the 
contact theorem of neutral polymers In our case the calculation of fix = 0) yields: 

pn.{w) = -^ 2 ((%(0) - \vicfiO) - 4) 2 +p$[e»<°> - 1] + 2c 6 [cosh(j/(0)) - 1] (14) 



where y(0) = /3e - 0(O) is the reduced electrostatic potential at the symmetry plane (x = 0). The above expression is 

b 



obtained by inserting the equilibrium (Boltzmann) distribution of the small ions c^- = exp( = p/3eV ) ) back into the 



free energy, eqs. 
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Equations [T^ and |l3| for the force are valid for the planar geometry (Fig. 1). In some experiments (l3) where 
the disjoining pressure of thin liquid films is measured, the two surfaces are indeed parallel to each other and H(w) 
is measured directly. However, most experiments ]5|-|l2]| use the surface force apparatus where the force Hr is 
measured between two cylindrical surfaces of radii R with a 90° tilt between their major axes (see Fig. 2a). At small 
distances w -C R the Derjaguin approximation fl30| relates the measured force to the excess free energy (as given by 
cq. [ll]) and not to its derivative (as given by eqs. p"2| , |l3| ) 

= 2nAF( W ) (15) 

R 

For clarity purposes we denote the force per unit area acting between two infinite flat surfaces as H(w) (eqs. (L2L[L3|), 
and the absolute force acting between two cross cylinders as ILr(w). 

The Derjaguin approximation can also be used to calculate the interaction between two spheres of radii R at small 
distances w <C R (Fig. 2b) as is the case for colloidal suspensions. The inter-colloidal force is related to the interaction 
free energy of two flat surfaces but with a different numerical factor: 

^ = kAF(w) (16) 



III. NUMERICAL RESULTS 
A. Low Salt Concentration 

Concentration Profiles: In order to solve the mean-field equations (eqs. ||J|), we use a minimal squares scheme 
in which the spacing between the two surfaces is divided into N ~ 100 intervals. An error functional which sums the 
squares of the local errors in eqs. is minimized with respect to the values of (f> and ip at the discrete grid points, 
under the constraint of the boundary conditions eqs. f?],|]. 

Typical solutions are presented in Fig. 3. The polymer is positively charged (p — 1) and attracted to non- adsorbing 
surfaces held at a constant negative potential (ip s < 0). Here we focus on the weakly screened limit (low salinity). The 
effect of screening can be estimated by comparing the inter-surface separation w with the Debye-Hiickel screening 
length kJ 1 defined by n 2 = QirlBCb where Ib = e 2 /eksT is the Bjerrum length equal to 7A for aqueous solutions at 
room temperature. At low salt concentration, kJ 1 3> w/2, screening is weak and plays only a minor role, whereas at 
high salt concentration, k^ 1 <C w/2, screening strongly reduces the Coulomb interactions in the adsorbed layer. In 
Fig. 3, the solution contains a small amount of monovalent salt (cb = ImM) so that the electrostatic screening length 
kJ 1 ~ 100A is larger than the inter-surface distance which varies between 7 and 40A. 

In Fig. 3a the reduced electrostatic potential y(x) = /3eip(x) is plotted as a function of the position x between 
the two surfaces. The reduced monomer concentration 4> 2 {x)/4) 2 is shown in Fig. 3b. Despite the fact that the 
surface potential is not very high, y s = —2.0 corresponding to %j) s ~ — 50mV, the adsorption is quite strong and 
the concentration in the gap between the two surfaces can increase by three orders of magnitude above the bulk 
concentration. The adsorption here is purely electrostatic since the only source of attraction is due to the electrostatic 
boundary conditions. A neutral polymer in similar conditions will not adsorb to the surfaces. 
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At small inter-surface distances the adsorbed polymers consist of a single layer extending from one surface to the 
other and the potential is negative everywhere in the gap. As the surfaces are drawn away from each other, first the 
amount of adsorbed polymer grows rapidly and then the adsorbed layer separates into two distinct layers near the 
two surfaces. In addition, the potential changes sign and becomes positive in the central region. By integrating the 
Poisson equation from the surface to the point where the potential changes sign, it can be easily shown that this sign 
reversal is due to an over-compensation of the surface charges by the layer of adsorbed monomers. 

The over-compensation is specific to charged polymers and does not appear in the Poisson-Boltzmann formalism 
for small ions (regular electrolytes). Physically, charged monomers which adsorb close to the surface are connected to 
other monomers which reside at some larger distance. In our model, the polymer chains resist fluctuations on length 
scales smaller than the Edwards correlation length £g ~ a/y/vc m of neutral polymers and thus over-compensate the 
surface charges and cause the potential to reverse sign. This over-compensation is much more pronounced when salt 
is added to the solution as will be discussed below. In the central region where the sign of the potential is opposite 
to that of the surface potential, the concentration of negative ions is larger than that of positive ions. We stress that 
when all contributions to the charge density are considered (polymer and small ions), the charges in the solution 
exactly balance the surface charges (see Fig. 5). 

At yet larger distances (w > 20A for the physical parameters of Fig. 3) the two adsorbed layers do not change any 
more. This occurs when the inter-surface distance w is larger than twice the width of the adsorbed layers. The two 
surfaces are almost decoupled and single surface adsorption is recovered. The polymer concentration between the 
two adsorbed layers is small and comparable to the bulk concentration. As long as the screening length kJ 1 is larger 
than the distance w, the electrostatic potential is nearly a constant (e.g. y ~ 0.2 in Fig. 3). At even larger distances, 
w > k~ , the effect of screening will show up, and the mid-plane potential will gradually decay to zero. 
Polyelectrolyte Adsorption: In Fig. 4a the total amount of monomers (per unit area) adsorbed between the 
two surfaces T(w) is plotted as a function of the inter-surface distance w for two charge fractions p = 1 (solid curve) 
and p — 0.2 (dashed curve). Similarly, in Fig. 4b the average reduced monomer concentration ((f> 2 /<j) 2 } is plotted as a 
function of w. Three regimes can be distinguished in accord with the findings presented on Figs. 3-4: at very short 
distances (w ~ 5 A) the confinement of the polymer to a narrow slit competes with the electrostatic attraction of the 
charged monomers to the surface and avoids strong adsorption. The polymer is not expelled totally from the gap 
between the two surfaces but its concentration is of the order of the bulk concentration. 

As the surfaces are taken further apart, the adsorption increases rapidly until it reaches its maximal value. At 
this point the average concentration can be three orders of magnitude higher than the bulk concentration. At larger 
distances the two surfaces decouple from each other and the adsorbed amount decreases towards a saturation value. 
At this stage the system can be described as two independent layers adsorbing onto the two surfaces. The saturation 
value of r is approximately twice the adsorbed amount to a single surface T(w — ► oo) = 2IV In a preceding work 
p6[ we investigated adsorption onto a single surface and have shown that in the low salt limit Ti ~ If-y/p- This 
behavior is a result of two competing interactions: (i) the electrostatic attraction of the charged monomers to the 
surface which is proportional to p; and (ii) the Coulomb repulsion between charged monomers in the adsorbed layer 
which is proportional to p 2 . For strong polyelectrolytes the latter interaction dominates and the adsorbed amount Ti 
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increases as the fractional charge p decreases. This scaling behavior is in accord with the saturated values of Fig. 4a. 

In Fig. 5 the charge densities per unit area are plotted as functions of the distance w for the same sets of values as in 
Fig. 4. Using the single surface results we verify that indeed in the saturated regime of large inter-surface separations 
a p (oo) ~ pTi ~ y/p. Another observation which can be made from Fig. 5 is that the two charge densities almost 
balance each other. In the low salt limit these charge densities must cancel each other since the amount of salt is too 
small to play any significant role in neutralizing the solution. 

Free Energies and Forces: In Fig. 6, inter-surface force profiles were calculated for Cb = 10~ 6 M, corresponding 
to k s _1 ~ 3000A, and for two values of the polymer charge fraction p. In Fig. 6a, the excess free energy per unit 
area 2irAF(w), eq. [H], which is the physical quantity measured in SFA experiments, is plotted as a function of the 
inter-surface distance w. In Fig. 6b the force per unit area II(k;), eq. [l^, acting between flat surfaces is plotted as a 
function of the distance. This force can be measured directly in disjoining pressure experiments of thin films |fL3| . 

The two surfaces strongly repel each other at very short distances and attract at short distances. Note that for 
small ions such attraction is not present within the Poisson-Boltzmann approximation. For polymer chains the lack of 
translational entropy as well as the large correlation length £e enhance the effective attraction between the surfaces at 
short separations. As the charge fraction is lowered, the attraction becomes weaker and the length scale of attraction 
increases. These two effects can be explained by simple arguments which are presented in section IV. Attractive 
interactions have been observed experimentally || and were attributed to "bridging" of chains between the surfaces 
— a mechanism which is not present in our approach. 

A secondary repulsion appears at large distances but is too weak to be shown on a linear scale in our plots. However, 
this secondary repulsion can be made quite pronounced, in particular when the polymer surface excess is fixed at a 
large value. This is further discussed in Sec. V. For small ions, the origin of the repulsion is entropic whereas here it 
is due to over-compensation of the surface charges. 

B. High Salt Concentration 

Concentration Profiles: At high salt concentration the screening length kJ 1 is smaller than the inter-surface 
separation w. The effect of screening on the concentration profiles is demonstrated in Figs. 7 and 8, where a strong 
polyelectrolyte (p = 1) is adsorbed from a solution containing large amounts of salt, Cf, = 1M for which the screening 
length is kJ 1 = 3A. In Fig. 7 the reduced electrostatic potential and the reduced monomer concentration are plotted 
as function of the position x for a range of inter-surface distances w. As a result of screening the attraction of 
charged monomers to the surface is reduced considerably as compared to the low salt limit and the total amount of 
adsorbed polymer is approximately half. Despite the weaker adsorption, the qualitative behavior is similar. At short 
distances (w < 15 A in Fig. 7b) a single adsorbed layer exists between the two surfaces. At intermediate distances 
(15A< w < 40A) this layer separates into two strongly interacting layers and at larger distances (w > 40A) the two 
layers decouple from each other and only weakly interact with each other. 

When the separation is large enough (w > 15 A in Fig. 7) the adsorption is strong enough so that the charged polymer 
over-compensates the surface charges. The signature of this effect is that the electrostatic potential changes sign in the 
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central region as seen in Fig. 7a. This sign reversal strongly affects the small ion concentration c ± (cc) as demonstrated 
in Fig. 8. Since the concentration of the small ions follows a Boltzmann distribution c ± (r) = exp(=F/3eV0 the 
concentration of small negative ions c~ in the central region exceeds that of the positive ions c + as can be seen in 
Fig. 8 for w = 40A (solid curve). Altogether, the net charge (per unit area) in the central region will be negative 
ensuring an overall charge neutrality between the two surfaces. Note, that at small separations (e.g. w — 10A in 
Fig. 8) the potential is negative everywhere and the concentration of positive ions exceeds that of the negative ions 
everywhere within the gap. 

Polyelectrolyte Adsorption: Despite the effect of screening, the adsorption of strongly charged polyelectrolytes 
is strong as can be seen in Fig. 7b. The reason is that screening has two competing effects: on one hand it reduces the 
attraction of charged monomers to the surface which is the driving force for adsorption. On the other hand screening 
also reduces the monomer-monomer Coulomb repulsion between adsorbed monomers, thus allowing for more charges 
to accumulate near the surface. Hence, despite the fact that the range of the electrostatic interaction is reduced 
considerably, the average polymer concentration near the surface can be high (as long as p is not too small). 

In contrast with the low salt regime, here the small ions play an important role in balancing the surface charges. 
In Fig. 9 the different contributions to the charge densities per unit area are plotted as function of the distance uu. 
The negative curve corresponds to the (induced) surface charge density on the two surfaces 2a s . The positive curves 
correspond to the total amount of polymer charges a p (per unit area) which have adsorbed between the two surfaces, 
and the total amount of charge carried by small ions <r + + a~ . At short distances, w < 6A, the polymer contribution 
is small and the main contribution to the charge density is that of the small ions. This distance can also be regarded 
as a lower cutoff to the continuum theory employed here since the monomer size we employed is of the same order of 
magnitude (a = 5A). 

When the surfaces are taken further apart, the polymer contribution and the net contribution of the small ions 
saturate to a constant value. This occurs when the two adsorbing surfaces decouple from each other, and two distinct 
adsorbed layers build up on each of the surfaces. Unlike the low salt case where the contribution of the small ions 
is negligible and the charged polymers dominate the charge density, in the high salt regime the contributions of the 
small ions and the polymer are comparable in magnitude. For example, in Fig. 9, the polymer contributes about one 
third of the charge density (per unit area) between the surfaces while the small ions contribute the other two thirds. 
Free Energies and Forces: Screening has a pronounced effect on the inter-surface interactions as can be seen 
in Fig. 10, where the inter-surface interactions are plotted as function of the distance for a strong polyelectrolyte 
(p = 1) at high salt concentration. The general behavior here is similar to the low salt case as is seen in Fig. 6. As 
the amount of salt increases, the adsorption is reduced because of the electrolyte screening and the attractive forces 
become substantially weaker. Similar effects were also observed experimentally in SFA experiments ]^JTl[|. 

Since the forces are related to the mid-plane values of ip and <ft ( ec l- 0) > it is of interest to study the mid-plane values 
of ip(0)- In the inset of Fig. 10 the mid-plane value of the reduced electrostatic potential y(0) — (3etp(0) is plotted as 
a function of the inter-surface distance for the same profiles that are used in the calculation of the forces. As can be 
seen also in Figs. 3 and 7, the mid-plane potential is negative at short distances, changes sign to become positive and 
finally decays to zero. The inter-surface interaction changes from repulsion to attraction at about the same distance 
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where the mid-plane electrostatic potential changes sign. This observation can be explained by examining the various 
contribution to the local free energy at the mid-plane, eqs. ^j4[ At the mid-plane the squared gradient terms in <j> and 
tjj vanish. In addition, for strong polyelectrolytes the excluded volume and chemical potential terms are very small. 
The force TL(w) — —f(x — 0) is dominated by two terms. 



The first term is the contribution of the charged monomers and changes sign when y(0) changes sign. The second is 
the (repulsive) osmotic pressure of the small ions and is proportional to j/ 2 (0) at small values of y(0). It is clear from 
the above that as long as the mid-plane potential is small enough, y(0) <C p</> 2 (0)/cf,, the pressure is governed by the 
first term and will change sign from attraction to repulsion when y(0) changes sign. It now follows that for a positively 
charged polymer, negative (positive) mid-plane potentials lead to repulsion (attraction) in agreement with Fig. 10. 
Since the potential changes sign when the adsorbed monomers over-compensate the surface charge, we conclude that 
this over-compensation is responsible for the reversal in the sign of the interaction. 

Weak polyelectrolytes (p <C 1) do not adsorb as much as strong polyelectrolytes. The attractive forces are much 
weaker and are easily overpowered by the double layer repulsion of the small ions. Such an example is presented in 
Fig. 11, where the inter-surface forces are calculated for a low charge fraction (p = 0.1) at two high salt concentrations 
Cb = 0.25M and c b = 1M. In contrast with the case of strong polyelectrolytes, here the forces are repulsive over the 
whole distance range and decay on a length scale of kJ 1 . 



The fundamental difficulty in studying polyelectrolytes is due to the competition between short range interactions 
such as the chain elasticity and excluded volume interactions, and the long range electrostatic interactions. In a 
previous work |2q] we have studied polyelectrolyte adsorption to a single surface by separating the two competing 
length scales: (i) The adsorption length D, which characterizes the width of the adsorbed layer and (ii) the electrostatic 
screening length kJ 1 = (Stt^q,)" 1 / 2 assuming that c;, 3> p4>1- The screening length depends on the salt concentration, 
while the adsorption length depends on both electrostatic and non-electrostatic properties. 

The two length scales can be separated in two limits: (i) the low salt regime D <C k" 1 and (ii) the high salt 
regime D ^S> kJ 1 . The difference between the two regimes being the range of the electrostatic interactions. The main 
assumption in this approach is that the polymer profile near a single flat surface can be written in the form 



where h(z) is a dimensionless function normalized to unity at its maximum and C sets the scale of polymer adsorption. 
The free energy can then be expressed in terms of D and C while the exact form of h(z) affects only the numerical 
prcfactors. Minimization of the free energy with respect to D and C gives the single surface adsorption length D\ 
and the concentration scale C\. 

When two surfaces interact with each other, the single surface profile is affected by the presence of the other surface. 
As a result the shape of the profile changes with the separation w as demonstrated in Figs. 3 and 7. For example, 



/m ~ -p0 2 (O)j/(O) + 2c b [cosh(y(0)) - 1 



(17) 



IV. SCALING REGIMES 




(18) 
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at short distances the profile varies monotonically between the surface and the half plane, while at larger distances it 
becomes non-monotonous, until finally the two surfaces decouple from each other, and the adsorption to each surface 
reduces to the single surface behavior. 

As in the single surface case, it is advantageous to separate the different length scales. First, we compare the single 
surface adsorption length D±, with the inter-surface separation w. At large separations w/2> D\ the surfaces interact 
weakly and the polyelectrolytes recover the single surface profiles. On the other hand, at short inter-surface distances 
w/2 <C Di, the gap is too small for the polyelectrolytes to follow the single surface profile. In this limit the relevant 
length scale (eq. [l8]) is just D = w/2, since w/2 serves as a lower cutoff for D. When w/2 increases so that w/2 ~ £q, 
the profile becomes more complex and our main assumption is no longer valid. Nevertheless, as demonstrated below 
this simplified picture reproduces the main features that characterize the inter-surface interactions in those short 
distances, w/2 < D\. 

Furthermore, the effect of screening can be taken into account by separating the screening length from the two 
other length scales. Two opposite limits are considered (i) The low salt regime kJ 1 ^> D\ and (ii) The high salt 
regime kJ 1 <C D\. 

A. Low Salt Regime: kJ 1 3> Di 

In the low salt regime the screening length is much larger than the width of the adsorbed layer and the effect of 
the small ions on the structure of the adsorbed layer can be neglected. This assumption amounts to neglecting the 
entropic contribution to the free energy /i on s(r) (eq. |^) and the electrostatic energies of the small ions in / e i(r) (eq. ||). 

1. Large Distances: w/2 S> Di 

At large distances w/2 D±, the two surfaces are only weakly coupled. The structure of the adsorbed layer near 
each of the two surfaces reduces to the single surface profile, and the "decorated" surfaces interact through a weak 
double layer interaction. In the limit of large distances and low salt conditions one needs to address the question of 
the relative size of w/2 and both lengths are large. 

The free energy of an isolated adsorbing surface can be approximated by [p6f 

/3F«(C, D) = a^C - a 2 p\y s \CD + 4tt/WC 2 j D 3 + ^(3 2 vC 2 D (19) 

The first term is the polymer elastic energy (or connectivity) term, the second term is the electrostatic interaction of 
the monomers with the surface, and the third term is the Coulomb repulsion between the adsorbed monomers. The 
electrostatic terms can be derived by integrating the interaction of every pair of charged layers at distances x and x' 
from the surface, with charge densities (per unit area) da = pe(j) 2 (x)dx and da' — pe<f> (x')dx' , respectively. Finally, 
the last term is the excluded volume term and will be neglected here since at low salt concentration its contribution 
is important only for extremely weakly charged polyelectrolytes. 

The coefficients ai,a 2 ,Pi and j3 2 are numerical prefactors, which depend on the exact shape of the dimensionless 
scaling function h(z). These coefficients can be explicitly calculated for a specific profile by integrating the Poisson 
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equation without taking into account the small ion contributions. For the simplest monotonous profile, namely 
a linear profile h\{z) = z for < z < 1 and h±(z) = for z > 1, we get a± — l,a^ — 1/3, j3\ = 1/14 and 
/?2 = 1/5. For a non-monotonous parabolic profile, hi(z) = 4z(l — z) for < z < 1 and fi2(z) = elsewhere, we 
get ai = 16/3, a2 = 8/15, Pi ~ 1/9 and /?2 — 2/5. Another profile which we consider is an intermediate profile of 
the form h^{z;rf) — A/r/ 2 z(rj — z) for < z < 1 and /i3(z) = elsewhere, where 77 is a parameter. This profile is 
non-monotonous and has a finite value at z = 1, which corresponds to the symmetry plane between the two surfaces. 
Furthermore, the special cases rj 3> 1 and 77 = 1 reduce to the simple linear and parabolic profiles h\(z) and /laO^), 
respectively. The parabolic profile hi{z) is a good choice for an isolated adsorbing surface in contact with a bulk of 
low concentration, whereas hi(z) describes better interacting surfaces at small separation. The third profile h^z;rj) 
can be regarded as intermediate between the other two. We stress that our scaling results do not depend on the 
specific shape of the profile h(z). Only the numerical prefactors will change. 

The single surface free energy (eq. |l9|) can be minimized with respect to both D and C along the same lines as was 
done in ref. ^6|. This yields a length scale D\ characterizing the adsorption onto a single surface 

d ^pW (20> 

and a concentration scale 

ci ~ j^J ( 2i ) 

In the low salt limit screening effects can be neglected as long as the screening length kJ 1 is much larger than the 
adsorption length D\. This condition limits the low salt regime to 

„ P\Vs\ f00 , 
c b < - — 2 (22) 

Inserting the above expressions back in the free energy gives the single surface free energy (up to numerical factors): 

n l/2\„ 15/2 

ftp™ ~ ~ P Vsl (23) 
ha 

At distances larger than the adsorbed layer x > D\ the amount of polyelectrolytes is small and comparable to its 
(low) bulk value. Since pip 2 Cb (even in the low salt limit), the interaction at large distances can be simplified. 
The system can be regarded as a solution containing electrolytes only (no polyelectrolytes) between two effective 
surfaces positioned at the edge of the adsorbed layers x = ±{w/2 — D\). The effective inter-surface distance is now 
w e ft — w — 2Di and each surface is kept at a (reduced) potential ud — P&ipD which is much smaller in magnitude 
than the original surface potential \i/>d\ <C |f/' s |. 

In the absence of polyelectrolytes in the effective gap, the electrostatic potential between to charged surfaces can 
be obtained by solving the Poisson-Boltzmann equation J|(J : 

V 2 i/;(r) = — c b sinh(/3e^) (24) 
e 

The above equation can be obtained from eq. ^| in the no polyelectrolyte limit. After the differential equation has 
been solved with the appropriate boundary conditions (namely, ip — tpo), the repulsive free energies and inter-surface 
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forces can be calculated. In particular, eq. |4|can be solved analytically in the weakly coupled regime, n s 1 <C w c s/2, 
yielding the following expressions for the repulsive inter-surface interactions: 

/3AF cl = tanh 2 (y D /4)e- Ksl0rff ~ 16c b Kj 1 y%e- K '' Wetl (25) 

/Hid = 64c 6 tanh 2 (j/ D /4)e- K8 ^« ~ 16c 6 ^e- K ^« (26) 

Other electrostatic regimes exist in which eq. |i] can be solved analytically. Those regimes lie beyond the scope 
of the present study because the polymer adsorbed layer reduces substantially the electric potential. Unfortunately, 
our model is too simple to give an accurate estimate of yo which is a local property. For this purpose more refined 
models are required. 

2. Short Distances; w/2 < Di < kJ 1 

At short separations w/2 -C D± <C kJ 1 the relevant length scale in eq. |l8| is w/2 instead of D. Consequently, 
and due to the planar symmetry of the system, the polymer contribution to the free energy can be written as 
F p (C,w) = 2F^ ] {C,D = w/2) where F { p 1] {C,D) is the single surface free energy (eq. |l^) . 

The free energy is minimized now only with respect to C leading to 

- \Vs\ w 2 -w 2 min 

C -J^ ^ ( 27 ) 

where w^ in = 2a\a 2 /2>a2P\y s \- The condition that C be positive, limits the validity of eq. |2^ to distances larger than 
a minimal distance w m i„, while at shorter distances C — and the polymers are depleted from the region between 
the surfaces. We estimate w m ; n ~ 0.21?! and so eq. ^ is valid in the range O.lDx < w/2 < Di |3(|. 

The total amount of monomers (per unit area) adsorbed between the two surfaces is directly related to C and C\ : 



1 mm 



pw/2 ( W < W n 

T(w) = / 4> 2 (x)dx = \ a 2 wC(w) w min < w < 2D X (28) 

•^-W2 I 2a 2 D 1 C 1 2D 1 < w 



Similarly, the average reduced monomer concentration is given by 

T(w) 



[,, min 



1 1 W 4>1 



w < w n 

a 2 C{w)/4>l w min < w < 2D X (29) 



2a 2 D 1 C 1 /w<i>l 2L>i < w 



The adsorption properties T(w) and (</> 2 /^) are plotted in Fig. 12 as function of the inter-surface distance w for 
the same physical values of Fig. 4. Our results agree qualitatively with the numerical results of Fig. 4 and reproduce 
the three different adsorption regimes. 

Inserting the above expression for C (eq. |2?|) back into the free energy, yields: 

(3F P * _]gl > 2 -^) 2 (30) 
Lb w 

for w m in < w < 2D\. At distances shorter than w m in, the polyelectrolyte is depleted from the gap, and the inter- 
surface interaction is dominated by electrostatic repulsion Bq]. 
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The inter-surface force I± p is readily obtained by differentiating the free energy F p with respect to w 

S((3F P ) _ \y s \ 2 (w 2 - w^Jjw 2 - bw^J 
p 5w l B w 6 

A quantitative comparison with the numerical results is shown in Fig. 13 where the physical parameters are the 
same as in Fig. 6. In Fig. 13 the polymer contribution to the excess free energy AF p (w) — F p (w) — 2Fp^ and to 
the inter-surface force n p are plotted as functions of the inter-surface separation w. The single surface free energy 
F^p is calculated from eq. ^3[ Only small separations w < 2D\ are shown in the figure. For w > 2D\ the shape of 
the inter-surface profile is more complex and we do not have scaling arguments relating the polymer profiles with the 
force. However, we expect the polymer contribution to be small. Comparing Figs. 13 with 6, we note that our scaling 
results (the characteristic length scale as well as the characteristic force) are in good agreement with the numerical 
results for several p values. 



B. High Salt Regime: k 3 1 < Di 



In the high salt regime the screening length is much smaller than the adsorption length D±. The Coulomb inter- 
actions between the charged monomers and the surface and between the monomers themselves decay exponentially 
with the Debye-Hiickel screening length kJ 1 . 

Our calculation is based on estimating the polymer contribution to the forces as mediated by the small ions. One 
should bear in mind that the contribution of the small ions to the forces is no longer negligible and can explain the 
discrepancy between the numerical (exact) and the scaling results. It is hard to get an analytical estimate to the 
small ions contribution because their concentration depends on the polymer profile via the electric potential. 



e Distance: «)/2> D\ 

The free energy (eq. |l^) can be generalized by introducing cut-off on the range of the electrostatic 

interactions (similar to what was done in Ref. J2(| for the single surface case): 

f3F p {C 1 D)=a 1 ^C-a 2 p\y s \CK: 1 +^I3 1 IbP 2 k~ 2 C 2 D + ^I3 2 vC 2 D (32) 

The electrostatic cut-off appears in two places. In the second term, only the first layers up to a distance kJ 1 from 
the surface interact with the surface charges. In the third term, each layer interacts only with its neighboring layers 
in the range of K~ ■ The numerical values of the prefactors ai,a2,/?i and (3 2 can in principle be different from the 
low salt values. However, since the prefactors are only used in Fig. 14 to demonstrate the qualitative behavior we will 
arbitrarily set their values to be the same as in the low salt regime. 
Minimizing eq. p2 with respect to D and C yields 



PWs 

and 



2 1/2 

Dl ~ W „ (33) 
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P 2 \y s \ 2 /M 2 
(3ip 2 /c b + f3 2 v 



which are now also functions of the salt concentration q, and k s — (S^so,) 1 / 2 ■ The condition that the screening 
length is much smaller than the adsorption length amounts to 

„ P\Vs\ rqc -S 

in agreement with the boundary of the low salt regime (eq. |3^) . 
The single surface free energy is now: 

PFU , - ^ffl' 1 '"'' (36) 

As in the low salt regime, the two adsorbed layers interact electrostatically. However, since the screening length 
is much shorter than the adsorption length in the high salt regime, this interaction decays quite rapidly. We note 
that the high salt regime can be further divided into two sub-regimes depending on the ratio of the two terms in the 
denominator of C\ and Fp as was discussed in Ref. [ p6| . 

2. Short Distances: to/2 < Di 



At short distances, the relevant length scale in the free energy eq. 32 is w/2 instead of D (as in the low salt regime). 
The free energy can be minimized with respect to C yielding 

c _ pIVsIkI 1 w-w min 

PlP 2 / Cb + @2V W 2 

where w mm — a\K s a 2 /3ct2p\y s \- As in the low salt case, C is positive only for w > w m - ln ~ 0.2D±. At smaller 
separations the polyelectrolytes are depleted from the gap and the inter-surface interaction is dominated by the 
electrostatic repulsion. We also note that the validity of eq. |3^ requires that w ^ ZkJ 1 . 
The polymer free energy of interaction is now 

flP _ p 2 \ys\ 2 K7 2 {w-w min ) 2 

PlP /Cb + P2V W d 

and the inter-surface force 



p 2 \y s \ 2 K s 2 (w - w min )(w - 3w min ) 
Pip 2 / c b + (3 2 v 



(3tt V \ys\ (*, \"J — w-anufyuj — uw mvo .) , . 

P V h -~ a ! —. — (3 y ) 



The qualitative behavior described by eqs. |3^ and ^ is similar to that of the low salt regime (eqs. [50|, |3l]). The 
typical behavior for a strong polyelectrolyte is depicted in Fig. 14, where the same physical parameters of Fig. 10 are 
used. We note that the quantitative agreement between the numerical (Fig. 10) and scaling (Fig. 14) results is not as 
good as in the low salt limit. For example, the value of the minimum in the free energy is about three times smaller 
in Fig. 14 as compared with Fig. 10. Also the variation of w at this minimum with the salt concentration is weaker in 
the numerical results. As discussed above, the main source of discrepancy between the numerical and scaling results 
is the omission of the small ion contribution in the latter. 
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C. Discussion 



To summarize our results we present in Fig. 15, a schematic diagram of the different adsorption regimes. The 
dashed lines mark the single surface adsorption regimes in terms of the charge fraction p and the salt concentration 
c&. Three adsorption regimes can be distinguished: 

• The low salt regime q, -C p\y s \/8TrlBO< 2 ■ 

• The first high salt (HS I) regime c& 3> p\y s \/8irlBa 2 with weak polyelectrolytes p 2 -C vc&. 

• The second high salt (HS II) regime Cb ^> p\y s \/8nlBa 2 with strong polyelectrolytes p 2 ^> vcb- 

The shaded area in Fig. 15 marks the region in parameter space where the polymer contribution to the inter-surface 
interaction is comparable to or larger than the pure electrostatic contribution. The shaded area includes the low salt 
regime, a large portion of the HS II (high salt/strong polyelectrolyte) regime and a small portion of the HS I (high 
salt/weak polyelectrolyte) regime. The exact crossover lines depend, of course, on the numerical coefficients which 
are not included in our approximations. Nevertheless, the qualitative picture can be deduced from the diagram. The 
different behaviors (as depicted previously) can be demonstrated with the help of this diagram. The filled circles in 
the low salt regime mark the graphs of Fig. 6 and are well within the shaded area. The filled squares on the right 
border of the diagram (at p = 1) correspond to the graphs of Fig. 10 representing strong polyelectrolytes in the high 
salt regime. At higher salt concentration the system is closer to the boundary of the shaded area and the polymer 
attraction is weaker. Finally, when the ionic strength is high enough the attractive contribution is too weak to be 
observed. Weak polyelectrolytes in the high salt regime belong to the top left side of the diagram outside of the 
shaded area. In this regime, the electrolyte dominates the inter-surface interactions which are purely repulsive, as is 
indeed the case for the force curves of Fig. 11. These curves are represented in Fig. 15 by filled triangles. 

In the following, we briefly summarize our findings in the different adsorption regimes. In Sec. VI the findings are 
compared with experimental works which are reported in the literature. 

Low Salt Regime: In the low salt regime, the Debye-Hiickel screening length is much larger than the width of 
the adsorbed layer. As a result the electrostatic interactions of the charged monomers with the surfaces and their 
interactions with other monomers are unscreened. This leads to strong adsorption as can be seen from the numerical 
results shown in Figs. 3-4 and from the scaling results (eqs. ^9|) shown in Fig. 12. In addition, since the bulk 
concentration of the salt and counter-ions, c& and pcj) 2 respectively, is small, the charge density in the solution between 
the two surfaces is mainly due to the charged monomers. This is demonstrated in Fig. 5, where it is shown that the 
surface charges are balanced by the charged monomers. 

At large distances the adsorbed polymer forms two distinct layers on the two surfaces. The amount of polymer 
adsorbed between the two surfaces saturates to a constant value (Fig. 4) which is approximately twice the single 
surface adsorbed amount. As discussed in a preceding work p6[ , the width of the single surface adsorbed layer D\ 
and the single surface adsorbed amount Ti both scale as p~ x / 2 . The dependence on p is due to the balance between the 
attraction of the monomers to the surface which is proportional to p and the monomer-monomer Coulomb repulsion 
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which is proportional to p 2 . The fact that the adsorbed amount decreases when the polymer charge increases reflects 
the energy barrier for bringing a large amount of charged monomers to the vicinity of the charged surface. 

The two layers start to overlap when the inter-surface distance is about twice the width of the single surface 
adsorbed layer w ~ 2D±. Below this distance the adsorbed amount slightly increases (Fig. 4) and the two surfaces 
strongly attract each other (Fig. 6). Our scaling approach recovers the increase in the adsorbed amount (eq. p8| and 
Fig. 12) and the attraction of the two surfaces (eqs. |30| , [3l| and Fig. 13). The magnitude of the polymer contribution 
to the interaction free energy scales as 

B l/2|y I 1/2 

PAF P ~ P Vsl (40) 
l B a 

This energy scale should be compared with the electrostatic interaction energy which scales as 

(3AF et ~ CbK~ 1 \y s \ 2 (41) 
The condition that AF p will be at least comparable to AF e i limits the salt concentration to 

l B a z 

in agreement with the boundary of the low salt regime (eq. |2^) . 

If the inter-surface distance is further reduced, the entropy loss due to the confinement of the polymer to a narrow 
slit pushes the polymer out of the gap between the two surfaces. This can be seen from the numerical results (Fig. 4) 
and also from eq. where w m { n is the minimal distance below which the polymer is compelled to leave the gap. 
High Salt Regime: In the high salt regime the Debye-Hiickel screening length kJ 1 is much smaller than the 
width of the adsorbed layer. As a result the range of the electrostatic interactions is much shorter and each charged 
monomers interacts only with monomers at a distance smaller than k~ . 

The limiting behavior at large distances depends strongly on the charge fraction p. For weak polyelectrolytes where 
p is small (regime HS I) the monomer-monomer Coulomb repulsion which is proportional to p 2 is negligible and the 

1 /2 

single surface adsorbed amount Ti scales as p/c b . On the other hand, for strong polyelectrolytes where p is large 
(regime HS II) the monomer-monomer Coulomb repulsion is dominant and the single surface adsorbed amount Ti 

1 /2 

scales as c b /p. The latter behavior is similar to that of the low salt regime with a different p dependence. At higher- 
salt concentration the adsorbed amount increases as the monomer-monomer Coulomb repulsion at the adsorbed layer 
is screened out. 

In the high salt case, the contribution of the small ions to the charge density can not be neglected. As seen in 
Fig. 9, about two thirds of the surface charge are balanced by small ions and only one third by charged monomers. 
Another interesting aspect of the interplay between charged polymer chains and small ions is the spatial distribution 
of charges between the two surfaces. The polyelectrolytes are strongly adsorbed on the surface resulting in a sign 
reversal of the potential: it is negative at the surface and becomes positive at a distance of 6 — 7A. In order that the 
system will be overall neutral the central region between the two surfaces has an excess of negative ions as seen in 
Fig. 8. At short distances (less than 12 — 14A) this effect disappears and the potential is negative everywhere in the 
gap. As seen from Fig. 10 and eq. ^ the sign reversal of the mid-plane potential y(x = 0) is accompanied by a sign 
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reversal in the inter-surface force H(w). The potential sign reversal can also be seen in the low salt regime. Since the 
concentration of small ions is very small, the over-compensation effect is weak. 

When the two surfaces are brought closer together, w < 2Di, these layers start to overlap (See Fig. 7) and the 
adsorbed amount slightly increases. At this separation, strong polyelectrolytes induce strong attraction between the 
two surfaces (Fig. 10). The polymer contribution to the attraction can be estimated from our scaling approach to be 

(3AF P ~ f^' 1 3 , s (43) 



Following the low-salt discussion we compare this interaction with the electrostatic interaction energy (eq. 41). For 
weak polyelectrolytes p 2 <C vc b (regime HS I) the polymer contribution dominates for 

cl < f¥ (44) 
while for strong polyelectrolytes p 2 3> vcj, the polymer contribution is dominant at low salt concentration 

c b « P M (45) 



At very short distances w < iu m m (eq. 37), the polymer is depleted from within the gap as can be also seen in 
Fig. 9. 

V. IRREVERSIBLE ADSORPTION 

So far in this study we have assumed that the adsorbed layer is in thermodynamic equilibrium with a bulk reservoir. 
Hence the total amount of adsorbed monomers can vary and is determined by the free energy minimization (eqs. [l]- 
However, in physical systems the energetic barrier for detaching an adsorbed chain from the surface can be much 
larger than the thermal energy ksT. As a result the relaxation times towards equilibrium can be much larger than 
the experimental time scales. In this case, one can consider the amount of adsorbed monomers as fixed. 

In this section we study the effect of irreversible adsorption by excluding the possibility of polymer exchange with 
the reservoir while keeping the total amount of monomers adsorbed between the two surfaces fixed. For simplicity we 
limit ourselves to the low-salt case. The results can be generalized to the high salt case. In principle, the state of the 
system is determined by minimizing the free energy (eqs. [l]-|]) under the constraint that 

r w/2 

</) 2 {x)dx =r (46) 

-w/2 

where Tq is the predetermined value of the amount of polymer between the surfaces. This constraint can be introduced 
through a Lagrange multiplier A replacing the chemical potential term in F so that the functional to be minimized 
becomes 

(3F = 0F - \(^J 1 (j) 2 {x)dx -r ^j (47) 

The self-consistent field equation now reads 

a 2 

yV 2 (Xr) = + P (f>(3e^ - \(j> (48) 
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while the modified Poisson-Boltzmann equation (eq. ||) is not affected by the irreversibility of the adsorption process, 
because the counter-ions are still free to exchange between the reservoir and the adsorbed layer. 

Equations || and ^8| are solved under the constraint of eq. ^ where A is adjusted to give the desired value of To- 
Typical force profiles calculated numerically are presented in Fig. 16a, where the free energy 2irAF is plotted as a 
function of the inter-surface distance w. For comparison we plot on the same graph the equilibrium free energy (solid 
curve) of Fig. 6. The dotted curve was calculated for the same physical values with the additional constraint that the 
total amount of monomers adsorbed between the surfaces is fixed to the equilibrium value of single surface adsorption 
(or equivalently two surfaces held at large distances). This value is defined as r sat . As can be seen from Fig. 4a for 

P = i, r sat = o.onA- 2 . 

The free energy difference between the solid and dotted curves is quite small and appears only at short distances 
when the adsorbed amount in true equilibrium starts to deviate from its saturated value (see also Fig. 4). Different 
values of Tq are shown on Fig. 16a. For low values of Tq (e.g. Tq = r sat /2) the attraction is weaker, while for higher 
values (e.g. Tq = 2r sat ) the attraction is much stronger. 

In order to understand better this behavior we return to the scaling approximation where we consider first the case 
of irreversible adsorption on a single surface. We assume that a fixed amount T = To/2 is adsorbed to the surface. 
Since the adsorbed amount T is related to the length scale D and the concentration scale C through r = (X2CD, it is 
possible to express C in terms of D and thus the free energy (eq. [lj]) in terms of D and Tq. Neglecting the excluded 
volume term, equation [l9] now becomes: 

0FW{r o ,D) = ^^r - l P \y s \To + ^l BP 2 T 2 Q D (49) 



Minimization of the free energy with respect to D yields 



1/3 



^"K&r.) (50) 

Substituting D\„ back in the free energy gives 

^W^-^l^lro+a^V 73 ^ 73 (51) 

where we have omitted some numerical coefficients. The first term in eq. |5l] is simply the interaction energy of the 
charged monomers with the surface while the second term is a balance between the monomer-monomer Coulomb 
repulsion and the chain elasticity term. 

When the two surfaces are interacting with each other, we can write an explicit expression for F p (w) at small 
distances (w -C Mr) by replacing the adsorption length scale D with w/2. The free energy then becomes: 

(3F p (w) = |^4 r ° " fllZ-Po + ^Ibp 2 T 2 w (52) 

The interaction free energy is now AF p = F p (w) — 2F p 1 \ The surface interaction term cancels out and we are left 
with three terms: two positive terms from cq. |52] and a negative term from eq. ^l|. The latter scales as Tq^ 3 and is 
responsible to the attraction at short distances where the third term of eq. [52] becomes small. In Fig. 16b we plot 
the interaction free energy as calculated from eqs. pILm for the same physical values as in Fig. 16a. Although the 
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numerical coefficients can not be obtained accurately from this approach the qualitative behavior is in accord with 
the numerical results. 

VI. COMPARISON WITH EXPERIMENTS 

The experimental studies with the Surface Force Apparatus p|-|l2fl focus mostly on the repulsive interactions between 
adsorbing surfaces. Due to the limitations of the experimental technique, the attractive interactions at short distances 
appear as jumps in force-distance profiles. Nevertheless, some of the qualitative features can be deduced from the 
experiments and agree with our findings. 

Luckham and Klein || have measured interactions between mica surfaces in the presence of poly-L-lysine which 
is a strong polyelectrolyte (p = 1) at two different salt concentrations Cb = ImM and Cb = 0.1M. The inter-surface 
forces Hr(w) (eq. |l5|) were measured at distances 50A < w < 1200A. The forces were always repulsive, decaying 
exponentially as function of w with decay lengths comparable to the Debye-Hiickel screening length k~ . These 
forces can be interpreted as the electrostatic repulsion between two adsorbed layers at distances larger than the width 
of a single layer (eq. ^5|). Significant deviations were found between the first approach where the two surfaces are 
brought close together and subsequent decompression-compressions cycles. The amplitude of the repulsion in the 
latter case was strongly reduced while the decay length of the force remained quite the same. In addition, the amount 
of polymer adsorbed between the two surfaces, as estimated by refractive index measurements, was much higher than 
in the initial measurements. Those effects demonstrate that the adsorbed layers are not always in equilibrium and 
that compression might lead to strong adsorption of charged polymers. The adsorbed amount remains high when 
the surfaces are separated from each other due to the high energetic barrier for desorption. The reduction in the 
inter-surface repulsion when the adsorbed amount is high is in accord with our numerical and analytical results for 
the case of irreversible adsorption (Sec. V). 

Marra and Hair M have adsorbed poly(2-vinylpyradine) (P2VP) between mica surfaces. The pH of the solution 
was such that the polymer was fully charged (p = 1) during the experiments. At low salt concentrations the forces 
were repulsive at large distances with an exponential decay which is consistent with the Debye-Hiickel screening 
length. Attraction of about — 7mN/m was detected at distances between 11 and 40A. In the presence of 0.01M 
NaCl the magnitude of the repulsive forces increased while the attraction was reduced to — 3mN/m and shifted to 
distances between 25 and 80A. At even higher ionic strengths (0.1M NaCl) the attraction disappeared altogether but 
an additional non-exponential contribution to the inter-surface repulsion was detected at distances between 60 to 
100A. These results confirm our findings that the effect of salt is to increase the adsorption length D\ and to reduce 
the polymer attractive contribution to the inter-surface interactions. 

Claesson and Ninham |Q] have adsorbed chitosan, a cationic biopolymer of glucosamine segments, between mica 
surfaces in the presence of 0.01 wt % acetic acid. The chitosan charge fraction was controlled through the pH of the 
solution 

1()-(p h -p k o) 

P = 1 + 10 -(pH-pK ) ( 53 ) 
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where Ko = 10 _pK ° ~ 10 -6 5 is the dissociation constant of the chitosan monomers. At low pH, the polymer is fully 
charged and the interactions were repulsive in the first compression when the surfaces were brought into contact. A 
double layer repulsion was detected at large distances (w > 100A) and strong steric repulsion at shorter distances. 
Upon separation, attraction was detected at distances around 20 — 25A. At pH=6.2 (p ~ 2/3) the repulsive double 
layer interaction disappeared altogether. The disappearance of the electrostatic double layer interaction indicates that 
the surface charges are exactly balanced by the adsorbed polymers. Attraction was detected at distances of about 
20 — 25A and strong steric repulsion at shorter separations. At pH=9.1 where the polymer was only weakly charged 
[p ~ 1/400) the double layer repulsion was again the dominant interaction. 

Dahlgren et al. ||] have studied the effect of salt concentration by adsorbing poly ((3-methacrylamido propyl) 
-trimcthyl-ammonium chloride) (MATPAC) between mica surfaces. Three salt concentrations were considered 
c& = O.lmM, 0.01M and 0.1M. The forces were repulsive at large distances and decayed exponentially with de- 
cay lengths of 170A, 30A and 11 A, respectively. The first decay length is smaller than the Debye-Huckel screening 
length of the corresponding salt concentration (kJ 1 ~ 300A) due to the contribution of the counter- ions. At higher 
salt concentrations this contribution is negligible and the decay lengths agree with the expected screening lengths. 
Attractive interactions were detected at short distances of a few nanometers. The magnitude of these attractive forces 
decreased as the amount of salt increased in agreement with the numerical results of Fig. 10. 

In another work, Dahlgren et al. || have studied the effect of both charge fraction and salt concentration on the 
inter-surface forces. Three different polyelectrolytes with different charge fractions p were used: MAPTAC (p = 1) 
and two copolymers AM-CMA-10 (p = 0.1) and AM-CMA-30 (p = 0.3) which were prepared using different ratios of 
neutral acrylamide (AM) segments and positively charged (2-acroyloxy ethyl) -trimethyl-ammonium chloride (CMA) 
segments. For a fixed charge fraction (p — 0.3) three different ionic strengths were compared: Cb — O.lmM, Cb = 0.01M 
and Cb — 0.1M. These experiments correspond to a vertical scan in Fig. 15. At the lowest ionic strength, the system 
is in the low salt regime and strong attraction is detected at intermediate distances 40 < w < lOOA. At the next ionic 
strength (c;, = 0.01M) the system is in the lower part of the high salt regime and weak attraction is still observed at 
distances below 60A. At the higher value of Cb = 0.1M no attractive interactions are observed and the electrostatic 
repulsions dominate the inter-surface forces. For a fixed ionic strength (q, = O.lmM) Dahlgren et al. have compared 
the inter-surface interactions for three different charge fractions: p = 0.1, p = 0.3 and p = 1. This set of experiments 
corresponds to a horizontal cut in Fig. 15. At the lowest charge fraction the repulsive interactions are dominant while 
for higher values of p attractions is observed at distances below w ~ 100A. 

Finally, the effect of ionic strength was studied separately by Dahlgren fll]| . Two polyelectrolytes, poly (2- 
proplyionyloxy ethyl) -trimethyl-ammonium chloride) (PCMA) and MAPTAC which have different molecular weight 
but the same charge fraction (p = 1) were studied. The monomers of these polymers are large and therefore the 
regime boundaries in Fig. 15 should be shifted to lower salt concentrations. Dahlgren has compared several types of 
multivalent salts at intermediate ionic strengths which are equivalent to Cb = 0.1,0.2,0.3 and 0.6mM. In the higher 
ionic strengths (0.3 and 0.6mM) no attractive interactions were observed. However, at lower ionic strengths attraction 
was observed at distances below w ~ 120A and w ~ 180A for Cb = O.lmM and Q, = 0.2mM, respectively. The ratio 
between the two lengths scales is approximately ~ 1.4 in agreement with the adsorption length scale in the high 
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salt regime D\ ~ c^ 2 (eq. [33|). Furthermore, the adhesion force was measured as Il a d ~ 170 mN/m and Il a d — 250 
mN/m in agreement with the force scale of strong polyelectrolytes in regime HS II F p ~ l/c\ /2 (eq||). 

VII. CONCLUSIONS 

In this paper we calculated numerically concentration profiles of polyelectrolytes between two charged surfaces and 
studied the inter-surface interactions as a function of the distance between the charged surfaces. Over-compensation 
of surface charges by adsorbed monomers was found to be strongly related to the reversal of inter-surface forces from 
repulsion to attraction at short distances where the two adsorbing layers strongly overlap. 

The effect of the polymer charge and ionic strength on the inter-surface interaction is studied by means of a 
simple variational approach. Three main regimes are found: (i) a low salt regime, q, <C pIi/sI/SttIbo 2 ', (ii) a high salt 
regime c;, ^S> p\y s \/8irlBa 2 for weak polyelectrolytes p 2 <C vcb (HS I); and (iii) a second high salt regime for strong 
polyelectrolytes p 2 3> vcb (HS II) . 

In the low salt regime, strong repulsion at very short distances is a result of the polymer depletion from the inter- 
surface gap. As the distance increases to w ~ a/p 1 ^ 2 , strong attraction is due to overlap of the adsorbed layers. 
Finally, when the inter-surface separation is larger than twice the adsorption length of a single surface, the two 
adsorbed layers separate and repel each other electrostatically. In the HS II regime the behavior is similar to the low 
salt one, with a modified length scale of interaction given by K s a 2 /p. On the other hand, in the HS I regime, the 
polyelectrolyte attractive contribution is too weak to generate a similar attraction at short distances. Consequently, 
the inter-surface interaction is repulsive with a decay length of kJ 1 . 

Some of the features described above are also present in the discrete lattice model of Bohmer et al. [[l8| . In particular, 
attractive interactions between equally charged surfaces were obtained numerically (Fig. 9 of Ref . [jl8| ) . This attraction 
was attributed to bridging by polymer chains from one surface to the other. The lattice model contains information 
regarding the fine details of the polymer chains which are absent in our model. On the other hand, the continuum 
approach is a convenient starting point for analytical approximations such as the scaling approach presented in this 
work. Attractive interactions were also obtained by Podgornik |27j] for the case of fixed adsorbed amount and without 
considering the nonlinear excluded volume interaction. For polyelectrolytes in a poor solvent, Chatellier and Joanny 
have obtained oscillations in the polymer concentration as well as in the inter-surface interactions. 
The model presented here takes into account the important Coulombic degrees of freedom within the frame work 
of the Poisson-Boltzmann formalism. We solve the coupled non-linear equations for the electrostatic potential and 
polymer concentration which allows consistent treatment of excluded volume effects as well as strong potential and 
surface charges (not the linearized Debye-Hiickel version). This allows us to consider cases where the Coulombic 
degrees of freedom are a major perturbation on the adsorption of neutral polymers. 

In the same time our approach has several limitations, some of which can be improved. The polymers chains are 
treated within a mean-field approximation which misses certain properties of polymer statistics such as the chain 
connectivity, stiffening of the persistence length, finite chain corrections and more specific conformations of polymers 
close to surfaces (loops, tails and trains). 
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On the other hand, the simple model we present offers a qualitative picture of polyelectrolyte chains between 
surfaces and suggests several scaling regimes. It can be further improved to take into account more realistic situations 
such as surface heterogeneities and other geometries, various charge distributions (quenched and annealed) on the 
chains ^4|,^5|,[37j , pH effects for acidic and basic polyelectrolytes |2(| as well as finite ion or monomer sizes [p8|— p0[ . 

In this work we have assumed constant surface potentials. One could also consider constant surface charges. In the 
low salt limit the behavior is expected to resemble the case of fixed amount of adsorbed polymer since the charged 
monomers are main source of charges which are able to neutralize the surface charges. In the presence of salt this is 
no longer the case and the behavior can differ considerable. 

It would be interesting to have thorough experimental results on the effect of the charge fraction p and salt 
concentration c& on the nature and magnitude of the forces and the corresponding length scales. We hope that our 
present work will encourage such systematic experimental studies. 
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FIGURE CAPTIONS 



Fig. 1: Schematic view of a polyelectrolyte solution between two parallel charged surfaces at a distance w from each 
other. The solution contains polyelectrolyte chains and small ions. The surfaces are kept at a constant potential. 

Fig. 2: (a) Two half cylinders with a 90° tilt between the axes, as used in the surface force apparatus to measure 
inter-surface forces. The radii R of the cylinders are of the order of 1 — 2 cm while w ranges down to a few Angstroms, 
(b) Two spheres of radii R at a distance w. Typically, colloidal suspensions contain particles whose radii are of a few 
microns down to hundreds of Angstroms, while the stability is determined by the balance of forces at much smaller 
distances. 

Fig. 3: Profiles of (a) the reduced electrostatic potential y = j3etp and (b) the reduced monomer concentration 
4> 2 1 4>\ as functions of the position x between the two surfaces. The profiles were obtained by solving numerically the 
differential equations (eqs. 0,0) for several inter-surface distances. For comparison, the different profiles are plotted 
on the same axis so that all mid-planes (x — 0) coincide. The surfaces are placed at different distances from the mid- 
plane as indicated by the filled squares. In the numerical examples in Figs. 3 to 11 we assume the following physical 
parameters: the polymer concentration is <j>1 = 10 _6 A -3 with an effective monomer length a = 5A and excluded 
volume parameter v = 50A 3 . It is immersed in an aqueous solution (e — 80) at room temperature (T = 300K) and 
the surfaces are kept at a constant potential y s — (3eip s = — 2. In addition, in the current figure the polymer charge 
fraction p = 1 and the salt concentration ct, = ImM. The different curves correspond to separations of w — 40A (solid 
curve); w — 20A (dots); w — 15A (short dashes) and w = 10A (long dashes). 

Fig. 4: Adsorption of polyelectrolytes in a low salt solution between two charged surfaces for two different polymer 
charge fractions, (a) Total amount of monomers adsorbed between the surfaces per unit area T and (b) the average 
reduced monomer concentration ((f) 2 / 4>b) as functions of the inter-surface distance w. The salt concentration is 
c b = 10~ 6 M. The different curves correspond to charge fractions of p — 1 (solid curve) and p = 0.2 (dashed curve). 

Fig. 5: Charge densities per unit area as function of the inter-surface distance w. The two positive (upper) curves 
correspond to the total amount of polymer charge adsorbed between the two surfaces per unit area a p = peT. The 
two negative (lower) curves correspond to the induced surface charge density on both surfaces, 2a s . The contribution 
of the small ions to the charge density is not displayed. The physical parameters and notations are the same as in 
Fig. 4. 

Fig. 6: Inter-surface interactions for polyelectrolytes between two surfaces held at a constant potential and in a low 
salt solution, (a) The excess free energy per unit area 2irAF as a function of the inter-surface distance w. The factor 
2ir is used in order to enable direct comparison with SFA measurements (see eq. [l5]). (b) The force per unit area II 
between the two surfaces as a function of the inter-surface distance w. The physical parameters and notations are 
the same as in Fig. 4. 



Fig. 7: Profiles of (a) the reduced electrostatic potential y = (3eip and (b) the reduced monomer concentration <fi 2 j§\ 
as function of the position x between the two surfaces. Same physical values and notations as in Fig. 3 except for a 
much higher value for the salt concentration, c\, — 1M. The different curves correspond to separations of w = 40A 
(solid curve); w — 20A (dots); w = 15A (short dashes); w = 10A (long dashes) and w = 7A (dots and long dashes). 

Fig. 8: Profiles of the small ion concentration c + and c~ as function of the position x between the two surfaces for 
two of the separations presented in Fig. 7: w = 40A (solid curve) and = 10A (long dashes). 

Fig. 9: Different contributions to the total charge density per unit area as function of the inter-surface distance w 
for highly charged polyelectrolytes (p = 1) in the high salt limit (c& = 1M). The different curves correspond to twice 
the induced surface charge density 2cr s , the total amount of polymer charge adsorbed between the two surfaces per 
unit area a p , and the total amount of charge carried by small ions (er + + cr~). 

Fig. 10: Intcr-surface interactions for highly charged polyelectrolytes (p = 1) at high salt concentration, (a) The 
excess free energy per unit area 2ttAF and (b) the force per unit area II between the two surfaces as function of the 
inter-surface distance w. The salt concentration is Cb = 0.1M (solid curve) and c;, = 1M (dashed curve). The inset 
shows the mid-plane values of the reduced electrostatic potential y(0) as a function of w. 

Fig. 11: Inter-surface interactions for weakly charged polyelectrolytes (p = 0.1) at high salt concentration, (a) The 
excess free energy per unit area 2ttAF and (b) the force per unit area II between the two surfaces as function of the 
inter-surface distance w. The salt concentration is q, = 1M (solid curve) and Cf, = 0.25M (dots). 

Fig. 12: Adsorption of polyelectrolytes in the low salt regime as calculated from eqs. |2^,|2{| (a) Total amount 
of monomers adsorbed between the surfaces per unit area T, and (b) the average reduced monomer concentration 
(4> 2 /4>t) as function of the inter-surface distance w. The physical values and notations are the same as in Fig. 4. The 
numerical prefactors of the intermediate profile h^{z, rj) with 77 = 3/2 were used in the calculation. The vertical lines 
denote the distance where w — 2D\. 

Fig. 13: Polyelectrolyte contribution to (a) the interaction free energy 2irAF p (eq. ^), and (b) the inter-surface 



force, lip (eq. 31), in the low salt regime as function of the inter-surface separation w. Same physical values, notations 
and units as in Fig. 6. The numerical prefactors of the intermediate profile ^3(2,77) with r\ = 3/2 were used in the 
calculation of the interaction free energy and forces, and the numerical prefactors of the parabolic profile were used 
in the calculation of the single surface free energies. 

Fig. 14: Polyelectrolyte contribution to (a) the interaction free energy, 2irAF p (eq. |38|), and (b) the inter-surface 
force, lip (eq. |39| ), in the high salt regime as function of the inter-surface separation w. Same physical values, notations 
and units as in Fig. 10. The numerical prefactors of the intermediate profile ha(z, rf) with 77 = 3/2 were used in the 
calculation of the interaction free energy and forces, and the numerical prefactors of the parabolic profile were used 
in the calculation of the single surface free energies. 



Fig. 15: Schematic diagram of the different regimes as function of the charge fraction p and the salt concentration 
Cf,. Three regimes can be distinguished: (i) the low salt regime D\ -c KjT 1 ; (ii) the high salt regime (HS I) D\ ^S> k^ 1 
for weak polyelectrolytes p <§C {cbv) 1 / 2 and (hi) the high salt regime (HS II) D\ ^> kJ 1 for strong polyelectrolytes 
p 3> (cfcf) 1 / 2 . The shaded area denotes the region where the polymer interactions are strong enough so that inter- 
surface attraction can be observed. The filled circles correspond to the numerical force-distance profiles of Fig. 6, 
filled squares correspond to the numerical profiles of Fig. 10 and the filled triangles to the profiles of Fig. 11. 

Fig. 16: The effect of irreversible polyelectrolyte adsorption on the interaction free energy. The excess free energy 
per unit area 2ttAF is plotted as a function of the inter-surface distance w. The graphs in (a) were obtained by 
solving numerically the mean field equations ||,^8| under the constraint of eq. The solid curve corresponds to the 
equilibrium interaction free energies and is the same as the solid curve of Fig. 6. The four curves were obtained for 
the same physical values as in Fig. 5 with p — 1 while the total amount of adsorbed monomers was kept at a constant 
value Tq. The different curves correspond to Tq — r sat = 0.01lA~ 2 (dots), Tq = 2r sat (small dashes), Tq = 3r sat 
(long dashes) and To = r sat /2 (dot-dash line). The graphs in (b) were calculated from the analytical expressions of 
the scaling approach (eqs. |5l|]5^ ). 
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